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Two  developments  have  combined  to  significantly  impact  the  growing  number  of 
applications  dependent  . pon  digital  terrain  elevation  data,  mathematical  terraj 
modeling,  and  minicomputer  growth.  Digital  representation  of  terrain  form  has 
previously  required  vast  amounts  of  mass  storage  with  the  relatively  slow  speeij 
data  access  assoc.iated  with  large  databases.  A technique  has  been  developed 
for  compact  digital  storage  of  eleimtton  data_wj;iich  also  decreases  the  dataj ^ 
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cessing speeds  and  digital  storage  capabilities  along  with  steadily 
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vation (iata  is  presented. 
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NEAR  REAL  TIME.  - APPLICATION  OF  DIGITAL  TERRAIN  DATA  IN 
■ MINICOMPUTER  ENVIRONMENT 


INTRODUCTION 

The  U.S,  Army  Engineer  Topographic  Laboratories  (USAETL)  is  a Corp 
of  Engineer  facility  whose  charter  Jncludes  research  ^nd  develop- 
ment in  the  areas  of  mapping,  charting;  remote  sensing  and 
geodesy.  USAETL  currently  employs  approximately  300  personnel  in 
its  four  main  research  and  development  organizational  elements: 
the  Topographic  Developments,  Computer  Sciences  and  Geographic 
Sciences  Laboratories,  and  the  Research  Institute. 

The  Automated  Cartography  Branch  (ACB)  is  a fourteen  person  element 
of  the  Topographic  Developments  Laboratory  and  is  tasked  with  the 
application  of  latest  computer  assisted  data  processing  techniques 
and  display  equipment  to  cartographic  problems.  The  ACB  conducts 
research  at  the  exploratory,  advanced  and  engineering  development 
levels.  Their  responsibilities  include  the  development  of  equip- 
ment and  techniques  for  automation  of  the  extraction,  color 
separation,  formatting,  storage,  retrieval  and  symbolization 
d-tsnlay  plotting  of  digital  cartographic  data  from  photographic 
and  other  source  material. 

A few  of  the  problems  considered  by  the  ACB  and  the  areas  of  concern 
of  this  paper  are  the  development  of  techniques  for  the  creation 
and  efficient  access  of  digital  cartographic  information  databases, 
and  their  effective  interface  with  user  application  software  routines. 

Reseanrh  by  the  ACB  over  the  past  few  years  has  resulted  in  a 
promi.:,..ag  answer  to  these  problems  for  terrain  elevation  data,  a 
mathematical  terrain  description  based  upon  polynomials.  This 
technique  describes  the  form  of  the  earth's  surface  using  a sparce 
matrix  of  simple  polynomials.  The  main  advantages  of  this  approach 
are  accurate  terrain  description  with  very  low  digital  storage 
requirements  and  substantially  reduced  applications  processing 
requirements . 

The  development  of  the  polynomial  terrain  mode]  has  occurred 
simultaneously  with  substantial  advances  in  the  minicomputer 
industry.  The  ACB  has  recently  been  tasked  with  investigating 
the  affect  of  these  two  developments,  (namely  the  polynomial 
terrain  model  and  the  new,  more  powerful,  low  cost  minicomputers) 
on  digital  terrain  elevation  data-related  applications.  This 
paper  presents  the  results  achieved  In  the  first  three  months  of 
a study  by  the  ACB  directed  to  building  and  evaluating  the 
utility  of  a polynomial  database  on  a minicomputer.  A brief 
survey  of  the  types  of  applications  Itnpliinented  to  date  Is  presented 
for  the  purpose  of  conveying  to  the  reader  a feeling  of  the  speed 
and  versatility  of  digital  terrain  elevation  related  applications 
that  are  readily  achievable  in  near-real  time  on  a minicomputer. 

An  overview  of  the  polynomial  terrain  model  approach  is  presented 
followed  by  a brief  description  of  the  test  data  and  facility.  The 
results  achieved  thus  far  are  then  Included  in  the  form  of  a series 
of  minicomputer  generated  graphics. 
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The  polynomial  matrix  format  consists  ol  a sparce  grid  of  locally 
determined  and  locally  valid  node  polynomials.  These  polynomials 
are  called  node  because  there  is  one  polynomial  centex*ed  over  each 
node  In  the  sparce  grii'.  The  processes  involved  In  their  generation 
and  use  (as  well  as  this  format’s  Important  characteristics)  will  be 
discussed  In  this  section  In  chronological  order. 

The  output  of  the  automatic  stereocompl latlon  equipment  Is  a computer 
magnetic  tape  containing  numerr.us  contiguous  DTM's,  each  In  sequential 
profile  format.  Batch  processing  software  Is  utilized  on  a large 
scale  digital  computer  to  mosaick,  transform,  and  regrld  these  DTM's 
into  a single  consistent  DTM  which  covers  the  area  of  interest  In 
the  desired  map  projection  or  cartesian  co-ordinate  system.  Further 
processing  may  be  utilized  to  edit  the  DTM  based  upon  separately 
digitized  planlmetrlc  detail,  such  as  lake  boundaries,  rivers, 
drainage,  and  rldgellnes.  It  is  at  this  point  that  the  polynomial 
modeling  process  begins.  Profile  subsets  of  the  DTM  are  sequentially 
defined  and  then  square  ground  areas  are  sequentially  defined  along 
each  profile  subset.  Each  square  ground  area  then  received  the 
follo^ring  processing.  (Note  that  since  each  subarea  receives  iden- 
tical, independent  processing  this  algorithm  contains  a high  degree 
of  parallelism  that  could  be  taken  advantage  of  on  one  of  the 
parallel  processor  digital  computers,  such  as  the  Goodyear  STARAN) . 
First,  the  subarea  Is  approximated  with  a "node"  polynomial  using 
the  least-squares  criteria.  At  this  point,  surface  slope  filtering 
and  sliqtle  statistical  measures  are  utilized  to  locate  extreme 
errors  In  the  data  and  to  measure  the  goodness  of  fit.  The  effects 
of  the  bad  data  on  the  approximation  are  efficiently  removed  through 
the  application  of  the  Inverse  of  the  Kalman  Filter  algorithm. 

Program  modifications  are  underway  so  that  the  length  of  the 
polynomial  can  be  automatically  varied  to  achieve  the  desired 
goodness  of  fit  and  error  distribution  bast'd  on  the  statistical 
analyses.  Independent  of  the  sl^e  of  the  square. area.  Currently, 
the  local  polynomials  are  the  most  efficient  four  coefficient 
polynomial , 


Z • Cy  + + 027  + C3xy 

the  simplest  non-trlvlal  case.  Th^  size  of  the  square  area  Is  chosen 
so  that  this  model  accurately  portrays  the  terrain.  Note  that  this 
order  node  polynomial  when  combined  with  the  weighting  functions, 
allows  for  the  modeling  of  terrain  forms  with  a wavelength  as  low  as 
the  node  polynomial  spacing.  When  these  processes  are  complete, 
the  node  polynomial's  coefficients  can  be  written  out  as  the  database. 

The  first  step  necessary  in  further  use  of  the  node  polynomial 
approximations  Is  their  combination  with  weighting  functions.  An 
algorithm  defines  the  final  model  over  each  square  area  defined  by 
the  centroids  of  four  neighboring  node  polynomial  approximations 
in  the  square  grid  as  a single,  different  simple  function  using  the 
weighting  functions  and  the  four  node  polynomials.  The  globally 
smooth  and  continuous  model  defined  results  in  a simple,  low  order 
polynomial  ratio  In  each  of  the  square  areas  defined  by  the  centroids 
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of  the  node  polynomials.  Grid  generation  and  analytic  solution  for 
the  contour  lines  arc  very  efficient  because  of  the  simplicity  of 
the  global  model. 


The  general  characteristics  and  strengths  of  this  digital  terrain 
format  can  be  summarized  as  follows:  (1)  automatic  model  accuracy 
definition  through  the  use  of  approximation  theorems,  statistical 
analyses,  and  variable  length  polynomials.  (2)  the  modeling 
algorithm  is  aroeanable  to  implementation  on  a parallel  Processor 
because  of  the  identical  Independent  processing  steps  fcr  the 
local  polynomial  subareas.  (3)  automatic  removal  of  large  data 
errors  via  surface  slope  filtering  and  inverse  Kalman  Filter 
algorithms.  (4)  efficient  model  format  requiring  minimal  core 
storage  because  of  the  sequential  and  independent  determination  of 
the  local  polynomials.  (5)  low  storage  overhead  since  the 
horizontal  position  of  the  local  polynomials  are  implied  by  their 
location  in  the  matrix.  (6)  a smooth  and  continuous  globally 
valid  model  through  the  use  of  WIT.  (7)  efficient  geographic 
access  of  the  database  because  access  to  a horizontal  position  is 
given  by  the  structure  of  the  matrix.  (8)  efficient  DTM  grid 
generation  due  to  the  local  low  order  polynomial  definition  of  the 
global  model. 


TEST  DATA  AND  FACILT  .Tf 

The  equipment  being  utilized  for  this  study  are  as  follows:  a DEC 
PDPll/45  minicomputer  with  128k  of  core  memory,  a 44  million  word 
disk  drive,  and  a Tektronix  4014  line  drawing  storage  CRT  terminal. 
The  test  area  is  the  600+  square  kilometers  covered  by  the  Cache, 
Oklahoma  1:50,000  mapsheet.  Normally  represented  by  over  2 million 
-16  bit  elevation  values,  this  area  was  stored  on  the  minicomputer 
disk  in  polynomial  format  using  only  12,000  polynomials  requiring 
only  32  bits  each — a compaction  ratio  of  over  80  to  1. 

Al.l  the  following  graphics  were  produced  on  the  minicomputer- 
controlled  CRT,  based  upon  interactive  user  commands.  Application 
routines  are  flexible  in  that  the  graohlcs  shown  can  be  produced 
at  anv  user  specified  position,  scale,  viewing  angle,  etc. 
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